TIME  COMPARISON  VIA  OTS-2 


G.de  Jong  and  R.Kaarls,  Van  Swinden  Laboratory (VSL) ,  Delft, 
Netherlands 

D.Kirchner  and  H.Ressler,  Technical  University (TUG) ,  Graz, 
Austria 


ABSTRACT 


Since  1980  time  comparisons  are  carried  out  via  OTS-2 
between  the  Technical  University  Graz  (Austria)  and  the 
Van  Swinden  Laboratory  Delft  (Netherlands) . 

The  method  has  been  based  on  the  use  of  the  synchronization 
pulse  in  the  TV-frame  of  the  daily  evening  broadcasting  of  a 
French  TV-prograrnme  to  Northern  Africa. 

Corrections,  as  consequence  of  enanqcs  in  the  position  of  the 
satellite  coordinates  are  applied  weekly  after  reception  of 
satellite  coordinates  from  EEOC. 

A  description  of  the  method  is  given  as  well  as  some  of  the 
particular  techniques  used  in  both  the  participating  labora¬ 
tories  . 

Preliminary  results  are  presented. 

Starting  January  1982  the  experiment  will  be  supported  by 
the  EEC  and  enlarged  by  participation  of  NPL ,  PTB ,  DFVLR , 

LPTF  and  IEN . 


INTRODUCTION 


In  order  to  improve  ( inter-) national  time  comparisons  over  long 
distances,  as  much  as  possible  independent  of  atmospheric  dis¬ 
turbances  and  other  obstacles  between  the  participating 
institutes,  the  application  of  several  types  of  satellites  and 
pulse  transmission  systems  is  investigated. 

Although  several  satellite  systems  have  proven  to  give  possibi¬ 
lities  for  excellent  time  comparisons,  so  far  most  of  the  systems 
are  limited  from  the  point  of  view  of  continuous  availability, 
continuous  accessibility,  on-site  installation  possibility, 
high  costs  and  so  on. 

Therefore,  an  investigation  is  carried  out  after  the  application 
of  always  available,  on-site  receivable,  easy  accessible,  low 
cost  satellite  systems. 

The  experiment  described  here,  makes  use  of  the  OTS-2  satellite. 
The  first  stage  is  based  on  the  application  of  a  one-way  method 
by  making  use  of  the  TV-synchronisation  method.  In  this  case 
corrections  must  be  applied  for  the  changes  in  the  position  of 
OTS-2 

In  a  later  stage  a  partly  two-way  method  will  be  considered. 
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MEASURING  PRINCIPLE 


The  time  comparison  is  based  on  the  measurement  of  the  time 
of  arrival  of  the  same  agreed  event  in  the  signal  from  a  geo¬ 
stationary  communications  satellite  at  both  sites  A  and  B 
(figure  1). 

The  signal  used  is  a  TV-broadcast  from  a  groundstation ,  which 
is  relayed  by  the  transponders  on  board  of  the  satellite. 

At  sites  A  and  B  receiving  only  stations  have  been  installed. 
Therefore,  the  method  is  a  one-way  method. 

The  time  difference  of  the  clocks  at  A  and  B  can  be  expressed 
as  : 


T  and  T  are  the  times  of  arrival  of  the  same  reference 
A  B 

pulse  as  determined  on  the  time  scales  in  A  and  B 

respectively ; 

t  '  and  4'b  '  are  the  propagation  times  from  the  satellite  to 
A  and  B  respectively; 

t  "  and  t  "  are  the  delay  times  of  the  signal  in  the  stations 
A  B  7 

A  and  B  respectively; 

L  and  7?  are  the  ranges  from  the  satellite  to  A  and  B 
A  B 

respectively; 

a  is  the  speed  of  light. 

t  "  and  t  "  are  constant  and  can  be  determined 
A  B 

t  '  and  t  '  are  depending  of  satellite  movements  of  the 
A  B 

geostationary  satellite. 

In  order  to  achieve  high  accuracy,  the  geographic  position 
of  the  receiving  stations  and  the  position  of  the  satellite 
must  therefore  be  known  with  high  accuracy. 

Due  to  the  relative  small  baseline  from  A  to  B  and  the  large 
ranges  R  and  to  the  satellite  a  geometric  conf iguration 

exists  in  which  the  first  order  influences  of  the  satellite 
movements  are  cancelled. 
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OTS-2  SATELLITE 

The  satellite  used  is  the  Orbital  Test.  Satellite  (OTS-2)  , 
launched  by  the  European  Space  Agency,  or:  May  11,  1978  and 
positioned  at  10  °E  in  the  geostationary  orbit.  The  expected 
lifetime  of  the  satellite  extends  to  1983. 

In  the  communications  field  several  experiments  arc  carried 
out  with  OTS-2. 

Some  details  of  the  satellite  ana  its,  payload  are  give  in  the 
figures  2  through  5. 

The  time  comparison  experiments  are  done  by  making  use  of  a 
TV-signal  sent  from  Prance  to  Northern  Africa  during  the  evening 
hours . 

The  received  carrier  frequency  is  1  1  .682  Gllz.  Tne  TV-signal  is 
FM- modulated  and  occupies  a  bandwidth  of  36  MHz.  The  signal  uses 
energy  dispersal  following  CCIR- recommendations . 

OTS-2  ORBIT  DETERMINATION 

The  satellite  position  is,  calculated  by  the  European  Space 
Opera t. ions  Centre  (ESOC)  in  Darmstadt  ~RG)  . 

The  ranging  data  are  determined  by  the  ground  stations 
[Ticino  (Italy)  and — partially —  Villairar.ca  (Spain). 

From  these  data  ESOC  calculates  the  OTS-2  orbit,  which  gives 
an  accuracy  of  about  2,6  km  (la)  for  the  posit. ion,  when 
manoeuvres  are  absent.  This  orbit  determination  is  carried 
out.  every  week  over  a  period  which  lasts  from  Wednesday  to 
Wednesday . 

Taking  into  account  the  geometry  of  the  position  of  TUG, 

VSL  and  OTS-2  the  next  table  shows  the  consequences  of  position 
uncertainty  of  OTS-2  on  the  accuracy  of  the  time  comparison. 

position  uncertainly  east-west  north-south  height. 

2  .  is  km  2 . 5  km  2  .  b  km 

resulting  time  error  17.9  ns  6b  ns  8  ns 

Table  1  :  Time  errors  resulting  from  the  uncertainty 

of  the  spacecraft  position 
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RECEIVING  STATIONS 
TUG,  GRAZ,  AUSTRIA 

The  timekeeping  station  of  the  Department  of  Communications 
and  Wave  Propagation  (INW)  of  the  TUG  is  located  at  the 
Observatory  Lustbuhel,  Graz. 

The  satellite  receiving  station,  which  is  also  managed  by 
the  INW  is  used  for  various  propagation  experiments  and  is 
therefore  rather  complex. 

In  figure  6  a  simplified  block  diagram  shows  the  configuration 
for  the  time  comparison. 

In  table  2  some  of  the  relevant  station  characteristics  are 
given. 

The  antenna  together  with  the  antenna  box  is  computer 
controlled  in  azimuth  and  elevation.  The  dish  can  be  heated 
and  the  antenna  box  with  the  RF-equipment  is  fully  aircon- 
ditioned. 

The  received  signal  is  converted  down  to  70  MHz  and  fed  into 
an  industrial  type  TV-receiver  with  a  bandwidth  of  36  MHz  and 
provided  with  a  threshold  extension  video  demodulator  and  a 
videoclamp , which  provides  more  than  40  dB  rejection  of  the 
energy  dispersal  wave  form. 

The  video  signal (IV  +  3  dB  into  75  ohm)  is  transmitted  by 
a  coaxial  cable  p^with  a  length  of  50  m  from  the  antenna 
building  to  the  time  keeping  station  in  the  observatory. 

The  station  clocks  are  HP  cesium  beam  frequency  standards. 

The  laboratory  is  further  equ.iped  with  several  Loran-C 
and  LF/VLF  receivers. 
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VSL ,  DELFT,  NETHERLANDS 

The  time-/ frequency  standards  laboratory  is  part  of  the 
Van  Swinden  Laboratory  of  the  National  Service  of  Metrology 
(Ministry  of  Economic  Affairs) . 

Figure  7  shows  a  simp .1.1. tied  block  diagram  of  the  receiver 
configuration . 

The  most  relevant  station  characteristics  are  given  in 
table  2. 

Tire  antenna  with  the  antenna  box  is  remote  controlled  in 
azimuth  and  elevation,  out  from  the  time-/ frequency  standards 
laboratory . 

The  antenna  is  placed  on  top  of  the  laboratory  and  is  connected 
to  the  time-/ frequency  standards  laboratory  by  a  coaxial 
cable  with  a  length  of  140  rn 


The  station  clocks  are  HP  cesium  beam  frequency  standards. 
The  laboratory  is  equiped  with  several  Loran-C  and 
LF/VLF  receivers. 
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THE  REFERENCE  POINT  IN  THE  SIGNAL 


The  common  reference  point/event  is  the  trailing  edge  of  the 

first  fieldsynchronization  pulse  of  the  video  signal. 

See  figure  8. 

Special  electronic  circuitry  has  been  developed  in  order: 

-  to  reject  the  DC-  offset  and  energy  dispersal  waveform; 

-  to  clamp  the  synchronization  pulses  to  zero  volt  level; 

-  to  measure  and  to  memorize  the  positive  white  peak  voltage 
of  the  video  signal ,- 

-  to  extract  the  reference  pulse  at  a  fixed  percentage  of  the 
maximum  positive  peak  voltage  of  the  video  signal; 

-  to  provide  the  positive  reference  pulse  with  an  amplitude  of 
more  than  2  V  into  50  ohm; 

-  to  suppress  the  reference  pulse,  when  the  video  signal  to 
noise  ratio  is  too  bad,  e.g.  when  there  is  no  TV-transmission; 

-  to  provide  a  positive  video-output  of  half  the  input  signal 
amplitude  into  75  ohm  for  a  video  monitor;  the  digital  sound- 
in-sync  pulses  are  removed,  so  a  stable  line  synchronization 
is  possible. 

The  delay  in  the  circuit  for  the  reference  pulse  is  about  10  ns 

and  is  very  stable;  the  pulse  to  pulse  jitter  being  only  de¬ 
pending  on  the  video  signal  to  noise  ratio. 
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MEASURING  SCHEME,  DATA  REDUCTION  AND  COMPUTATION 


Since  September  2,  1980  the  time  scales  UTC  (TUG)  and  UTC 
(VSL)  are  compared  by  making  use  of  the  TV-signals  transmitted 
by  the  OTS-2  satellite. 

Every  day,  beginning  at  lOhlOmin  00s  UTC  the  first  measurement 
is  made,  followed  by  one  measurement  every  two  seconds, 
up  to  a  total  of  20  measurements  (last  measurement  18hl0min38s) . 
This  series  of  twenty  measurements  is  repeated  every  hour  with 
the  last  series  ending  at  22h30min38s  UTC. 

Each  series  of  measurements  is  followed  by  a  Loran-C  measurement 
of  the  transmitter  Syit  (7970W) ,  which  is  one  of  the  slaves  of 
the  Norwegian  Sea  chain.  For  the  rest,  the  actual  number  of 
measurements  per  day  depends  mainly  on  the  presence  and 
duration  of  the  TV- transmission. 

The  reference  pulses  have  a  repetition  rate  of  50  Hz. 

The  time  difference  between  the  local  time  at  each  station 
and  the  reference  pulse  is  measured  with  an  accuracy  of  1  ns. 

The  measuring  data  are  exchanged  once  a  week. 

The  data  are  processed  in  the  following  way.  From  each  series 
of  20  measurements  a  mean  value  is  computed  by  an  iterative 
process  omitting  all  data  outside  an  originally  arbitrary  limit 
(+0.5  ps)  from  the  resulting  mean  value.  Moreover,  this  mean 
value  is  only  used  for  further  computation  if  the  number  of 
measurements  used  to  form  this  mean  is  greater  than  50  %  of  the 
available  measurements  (on  an  average  91  %  of  the  available 
measurements  are  used) .  The  1  a  standard  deviations  are  also 
computed  and  their  mean  value  is  about  30  ns . 

Also  the  1  0  standard  deviations  for  the  readings  of  each 
station  are  evaluated  after  removing  the  systems tics  by  a 
linear  regression  through  the  series  of  data  points  and  result 
on  an  average  in  a  figure  smaller  then  10  ns  for  the  TUG  station 
and  30  ns  for  the  VSL  station. 

Next.,  corrections  as  consequences  of  the  ranging  data  are 
applied.  This  leads  to  UTC  (TUG)  minus  UTC  (VSL). 

Only,  still  a  certain  offset  is  present  as  a  result,  of 
so  far  unknown  equipment  delays. 

The  determination  of  these  delays  has  still  to  be  carried 
out . 

To  identify  the  outliers  after  correcting  for  the  ranging  data, 
a  linear  best  fit  of  certain  length  (50  data  points,  i.e.  50 
hourly  mean  values  or  approximately  20  days)  is  applied. 
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If  the  data  point  in  question  deviates  by  no  more  than 
a  given  limit  {  +  0.5  jjs)  from  the  estimate  for  this 
point /which  is  given  by  the  above  mentioned  linear  regression 
through  the  last  50  preceeding  and  used  data  points, then  this 
point  is  used  and  becomes  the  last  point  in  the  next  calculated 
linear  best  fit. 

A  further  condition  is  that  only  such  data  points  are  used 
which  are  the  mean  of  at  least  10  single  measurements. 

For  further  processing  of  the  data  from  the  hourly  mean 
values  daily  mean  values  were  calculated  and  also  the  1  a 
standard  deviation  for  all  measurements  of  one  day 
(the  1  o  values  are  on  an  average  smaller  than  35  ns) . 

The  measurements  were  then  corrected  for  a  step  (1*6  x  lO-1^  ) 
in  the  relative  frequency  difference  of  the  both  time  scales. 
This  step  is  evident  from  the  OTS  and  LORAN-C  data. 

Thereafter  a  linear  best  fit  was  subtracted  from  the  data. 

The  associated  LORAN-C  values  were  processed  in  the  same  way. 

RESULTS 

Figure  9  shows  a  plot  of  the  hourly  mean  values  of  the 
readings  (TUG)  minus  readings  (VSL) . 

A  time  step  of  9  ys  in  the  UTC(TUG)  scale  at  MJD  44744 
(May  20,  1981)  is  removed  and  furthermore  the  values  for  UTC 
(TUG)  before  MJD  44564  (November  21,  1980)  were  replaced  by 
the  corresponding  values  of  a  second  clock (HP  5061 A  with  High 
Performance  Tube)  which  was  used  later on  to  generate  UTC (TUG). 
The  overall  slope  of  the  curve  is  caused  by  the  frequency  dif¬ 
ference  (about  7xlO“13)  of  the  involved  time  scales  UTC (TUG) 
and  UTC (VSL)  respectively.  The  characteristic  shape  is  due  to 
the  change  of  propagation  time  resulting  from  the  movements 
of  the  satellite. 

Figure  10  is  a  plot  of  all  points  which  meet  the  specifications 
in  order  to  be  taken  into  account  for  the  calculation  of 
UTC (TUG)  minus  UTC (VSL) . 

Also  plotted  in  figure  1-0  are  the  values  of  portable  clock 
visits  by  the  U.S.  Naval  Observatory  to  both  stations. 
Unfortunately  however,  the  clock  visits  to  both  the  stations 
were  not  carried  out  on  the  same  trip.  Therefore,  the  plotted 
values  are  extrapolations  using  the  clock  rates  to  obtain 
UTC (TUG)  minus  UTC (VSL)  from  UTC(USNO)  minus  UTC (TUG)  and 
UTC(USNO)  minus  UTC (VSL) . 


Fig.  11  shows  the  residuals  of  the  OTS-2  and  the  LORAN-C  data 
plus  arbitrary  constants.  The  third  trace  of  data  points  in 
Fig.  11  presents  the  difference  of  the  OTS-2  and  LORAN-C  data 
for  mean  values  of  10  days  (a  constant  is  added  to  the 
difference  so  that  the  first  value  is  equal  to  zero) , 

Also  indicated  are  the  differences  between  the  OTS-2  results 
(regression  line)  and  that  obtained  by  the  portable  clock 
visits  of  the  USNO  (  100  ns  and  -  43  ns  respectively) . 

No  attempt  was  made  to  correct  for  the  difference  in  the 
delays  of  the  receiving  equipments. 

The  accuracy  and  precision  of  this  result  is  mainly 
determined  by  the  knowledge  of  the  position  of  the  spacecraft. 
Figs.  12  and  13  show  a  portion  of  Figs.  S  and  10  on  an 
expanded  scale  to  demonstrate  this  effect. 

The  data  in  Fig.  13  are  corrected  for  the  frequency  difference 
of  the  two  time  scales  and  clearly  show  the  influence  of  the 
uncertainties  caused  by  the  ranging  data. 

Fig.  14  is  a  stability  plot  (two  sample  Allan  variance 
0^  (2,  x)  )  of  the  UTC(TUG)  -  UTC(VSL)  time  scales  compared 
by  TV  via  OTS-2  and  LORAN-C.  For  sampling  times  longer  than 
5  days  the  stability  of  the  OTS-2  TV-link  is  better  than  that 
of  the  LORAN-C  signal.  Also  plotted  in  that  figure  are  the 
a  (2,  t)  values  if  one  only  uses  the  data  period  from  MJD 
44646  to  MJD  44665  (see  also  Fig.  13) .  These  figures  give 
an  estimate  of  what  could  be  gained  by  time  comparisons  using 
the  TV-signal  transmitted  by  the  OTS-2  satellite  if  better 
tracking  data  would  be  available. 

CONCLUSIONS 

-  The  given  results  are  not  corrected  for  equipment  delays, 
so  the  data  have  a  certain  offset. 

Also  a  direct  transportable  clock  comparison  should  be 
carried  out  in  the  near  future. 

-  The  ranging  data  could  be  improved  by  the  application  of 
another  best  fit  method  by  ESOC. 

This  was  discussed  with  ESOC  and  will  be  carried  out  by 
them. 

-  From  one  orbit  determination  period  to  the  next  an  offset 
of  up  to  200  ns  can  be  seen,  especially  when  there  was  a 
spacecraft  manoeuvre  in  between  the  periods  of  computation 
(see  Fig.  13  :  MJD  44645) . 

-  The  receiver  at  the  VSL  could  be  improved  in  order  to 
achieve  less  scatter  in  the  readings. 

This  improvement  has  been  carried  out  recently. 
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-  The  Loran-C  data  seem  to  have  a  strong  seasonal  dependence. 

The  difference  between  the  Loran-C  and  OTS-2  results 
reached  700  to  900  ns  for  a  differential  path  length  over 
land  of  600  km. 

The  la  standard  deviation  for  the  period  of  one  year  is 
107  ns  for  the  OTS-2  data  and  270  ns  for  Loran-C. 

-  In  January  1982  the  OTS-2  experiment  will  be  extended  by 
the  participation  of  NPL,  PTB,  DFVLR,  LPTF  and  IEN. 

The  experiment  will  be  supported  by  the  EEC. 
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table  2 

Main  characteristics  of  the  receiving  stations 
TUG  VSL 

geographical  position  47°04'02.88"  N  52°00'01.85"  N 

1  5°29 ' 4 1 . 80"E  4°22  1  56 . 18"  E 

405.32  m  above  20.16  in 

sea  level 


antenna  diameter 

3  m 

3  m 

pre-ampli fier 

uncooled  parametric 
low  noise  amplifier 

Gci As  Fet. 

low  noise  amp li fie 

noise  temper- 

ature 

190  K 

317  K 

down  conversion 

two  stage  single  stage 

first  IF  (750j_  250)  MHz  (135  +  18)  MHz 
second  70  MHz 

antenna  gain 

48  dB  at  11  GHz 

47.5  dB 

system  noise  temper- 

ature 

370  K 

832  K 

figure-of -merit 

22,3  dB/K 

18,3  dB/K 

carrier  to  noise  ratio 

16,7  dB 

12  dB 

TV-receiver 

double  conversion 

single  conversion 

f:i  Iter  bandwidth 

36  MHz 

40  MHz 

time-  interval  counter 

2  ns 

1  ns  acc. 

data  storage 

magnetic  tape 

diskette 

The  geographical  coordinates  of  the  antennae  at  both  the 
stations  are  determined  by  Doppler  measurements  carried 
out  in  Western  Europe  during  the  last,  years.  The  calculations  were 
made  by  IFAG,  Frankfurt,  a. Main,  (FRG)  . 
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Overall  Characteristics  of  Communications  Payload  ( Beginning  of  Life) 
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Antenna  beam- Polarisation  Euro  A  — X  Euro  A  — Y  Euro  B  —  R  Euro  B  — L 

EIRP  at  beam  centre  19.5  19.5  t  31  2  31.1  dBW 
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QUESTIONS  AND  ANSWERS 


MR.  DAVID  ALLAN,  The  National  Bureau  of  Standards 

As  you  move  towards  10  nanosecond  level  of  accuracy,  I  noticed  in 
your  original  equations  the  sagnac  was  not  included,  and  that  be¬ 
tween  those  two  sites,  it  calculates  to  be  about  25  nanoseconds 
interestingly? 

DR.  KAARLS: 

What  was  the  correction.  I'm  sorry? 

MR.  ALLAN: 

The  sagnac  correction. 

DR.  KAARLS: 

Okay.  No,  it  was  not  included,  you  are  right. 

MR.  ALLAN: 

Yes.  And  especially  as  you  go  to  North/South  stations  where  the 
sagnac  will  be  very  small  with  respect  to  East/West  stations,  you 
will  get  incontinuity  or  discontinuities  if  you  don't  apply  that 
correction  as  you  get  down  to  the  10  nanosecond  region? 

DR.  KAARLS: 

Yes.  It  might  be  a  factor  of  importance,  yes.  Thank  you. 


369 


